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Abstract 

A set of exactly solvable one-dimensional quantum mechanical potentials is de- 
scribed. It is denned by a finite-difference-differential equation generating in the lim- 
iting cases the Rosen-Morse, harmonic, and Poschl- Teller potentials. General solution 
includes Shabat's infinite number soliton system and leads to raising and lowering op- 
erators satisfying q-deformed harmonic oscillator algebra. In the latter case energy 
spectrum is purely exponential and physical states form a reducible representation of 
the quantum conformal algebra su q (l, 1). 
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Lie algebras are among the cornerstones of modern physics. They have enormous number 
of applications in quantum mechanics and, in particular, put an order in classification of 
exactly solvable potentials. "Quantized", or g-deformed, Lie algebras (also loosely called 
quantum groups) are now well established objects in mathematics Their applications 
were found in two-dimensional integrable models and systems on lattices. However, despite 
of much effort quantum algebras do not yet penetrate into physics on a large scale. In 
this paper we add to this field and show that a g-deformed harmonic oscillator algebra @ 
may have straightforward meaning as the spectrum generating algebra of the specific one- 
dimensional potential with exponential spectrum. This result shows that group-theoretical 
content of exactly solvable models is not bounded by the standard Lie theory. 

Recently Shabat analyzed an infinite chain of reflectionless potentials and constructed an 
infinite number soliton system ||. The limiting potential decreased slowly at space infinities 
and obeyed peculiar self-similar behavior. We will present corresponding results in slightly 
different notations. We denote space variable by x and introduce N superpotentials W n (x) 
satisfying the following set of second order differential equations 

{W' n + W' n+ i + W 2 - W 2 +1 )' = 0, n = 0,...,N-l (1) 

where primes denote derivative w.r.t. x. Taking first integrals 

W' n + WUi + W 2 n - W 2 +1 = k n+1 , (2) 

where k n are some constants, we define N + 1 Hamiltonians 

2H n = p 2 + U n (x), p=-id/dx, (3) 

U (x) = W 2 -Wq + ko, U n+1 (x) = U n {x) + 2W' n {x). 

An arbitrary energy shift parameter fc enters all potentials U n (x). 

Notorious supersymmetric Hamiltonians are obtained by unification of any two successive 
pairs H n ,H n+ i in a diagonal 2x2 matrix 0. Analogous construction for the whole chain 
(0) was called an order iV parasupersymmetric quantum mechanics || ||. In the latter 
case relations ([!]) naturally arise as the diagonality conditions of the general (N + 1) x 
(N + 1) -dimensional parasupersymmetric Hamiltonian. We do not use here these algebraic 
constructions and consider operators H n on their own ground. 

If W n (xys do not have severe singularities then the spectra of operators (H) may differ 
only by a finite number of lowest levels. Under the additional condition that the functions 

^\x) = e-r^y)dv (4) 
belong to the Hilbert space £2 one obtains first N eigenvalues of the Hamiltonian Hq 

n 

H ^\x) = E n ^\x), £„ = §5>, n = 0,l,...,iV-l, (5) 

i=0 
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where subscript n numerates levels from below. In this case @) represents ground state 
wave function of H n from which one can determine lowest excited states of Hj,j < n, e.g., 
eigenfunctions of H are given by 

#0) oc (p + iW Q )(p + iWt) . . . (p + iW n _i) 4 n) . (6) 

Any exactly solvable discrete spectrum problem can be represented in the form 
Sometimes it is easier to solve Schrodinger equation by direct construction of the chain 
of associated Hamiltonians @ - this is the essence of so called factorization method [7- 
9]. For the problems with only N bound states there does not exist Wn(x) making V'o 
normalizable. For example, if Wn(x) = 0, then H n has exactly N — n levels, the potential 
U n (x) is reflectionless and corresponds to (N — n)-soliton solution of the KdV-equation. 

Let us consider potentials which support infinite number of bound states, N = oo. In this 
case one can derive from (|2|) the following differential equations involving only one derivative 
and a tail of W n 's 

oo 

W((x) + Wf{x) + ]T(-iy(2Wf +J (x) + k i+j ) = 0. (7) 

3=1 

A question of convergence of the infinite sum is delicate and requires special consideration 
in each case. Evident condition W 00 (x) = W^(x) = 0, which is still related to the soliton 
dynamics, is necessary for rigorous justification of (0). Here we shall operate with formal 
series and assume that initial chain (0) always may be recovered by adding (|7|) for i = n and 
i — n + 1. In order to find infinite number of superpotentials {W{\ from (^) one has to relate 
them to one unknown function via some simple rule. Following Ref. we take the Ansatz 



Wi(x) = q'Wtfx), (8) 



which yields the equation 



W'(x) + W\x) - 7 2 + 2j2(-l) j Q 2j W 2 (q j x) = 0, (9) 

j=i 

where 7 2 = — Y^=i{~^Ykj. Note that reality of superpotentials does not necessarily restrict 
parameter q to be real - this will appear later. From (§) it is easy to derive eqs. (0) and (fj) 
with 

ki+i = 7 2 (l + gV\ *>0. (10) 

The following computation 

oo oo 

7 2 = - E(-!) J % = 7 2 (1 + q 2 ) E(-l) i 9 2J = 7 2 (H) 

3=1 3=0 

shows that 7 2 is completely arbitrary parameter (an energy scale) and ( [TD[ ) is a self-consistent 
definition of the constants fcj. Derivation (|TT| ) is valid only at \q\ < 1, which was the 
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restriction of Ref.||, but if fl8|) and fllPf ) are taken as the basic substitutes for (§) then by 
definition j 2 is arbitrary and there are no essential restrictions on q up to now. 

Eq. (H) has certain relation to quantum algebras JTJ and corresponding g-analysis [|10|1 . 
In order to see this we first introduce a scaling operator T q obeying group law 

T q f(x)=f(qx), T q T p = T qp , T~ l = T q -i, T\ = 1. (12) 

Then (|]) can be rewritten as 

oo 

W'(x)-W 2 (x) = 7 2 -2]T(-l)VT g )W 2 (x) 

j=0 

= 1 2 -2{l + q 2 T q )- l W\x). (13) 

Multiplying ( |T3"[ ) from the l.h.s. by (1 + q 2 T q ) we obtain finite-difference-differential equation 
defining W(x) 

W'(x) + W 2 (x) + qW'(qx) - q 2 W 2 (qx) = 7 2 (1 + q 2 ), (14) 

which is nothing else than the first iteration of superpotentials. The whole infinite chain (Q) 
is thus generated by (|T4]). This observation removes ambiguities arising in ([J) due to the 
convergence problems. 

Let us try to find quantum mechanical spectrum generated by the self-similar potential 
U (x) associated to (|i~4]). Suppose that eigenfunctions (f|) are normalizable. Then potential 



U i+ i(x) contains one eigenvalue less than Ui(x), i.e. there should be the following ordering 
of levels 

n 1 -I- a 2 

E Q <E l <...<E 00 , E n = \Y,h = -|7 2 T^g 2n , (15) 

U 2 l-q 2 

where we chose undefined constant ko to be ko = — 7 2 (1 + (? 2 )/(l — q 2 )- At negative 7 2 it is 
not possible to fulfil the ordering and at positive 7 2 the parameter q should be real and lie in 
one of the regions \q\ < 1 or \q\ > 1. Taking the normalization 7 2 = u 2 \l — q 2 \/ (1 + q 2 ) and 
denoting \q\ = exp(±r//2), rj > 0, we arrive at exponentially small or large bound energy 
spectrum 

E n = T \u: 2 e^. (16) 
What type of potentials these spectra would correspond to? In order to know this one 



should solve equation (14). Then everything crucially depends on the normalizability of i/j^ 



in (f|) because all other wave functions ip^ are related to it by scaling. Normalizability is 
insured if W(x) is a continuous function positive at x — > +oo and negative at x — > — oo. 
Under such conditions W(x) has at least one zero and we choose corresponding point to be 
x = 0, i.e. W(0) = 0. Eq. ( |T4"D now automatically leads to W(—x) = —W(x) and below 
we restrict q to be semipositive. Let us find solution of ([14]) in the Taylor series form near 
the zero. Substituting an expansion W(x) = Y^Lx CiX 2% ~ x into (f|) we obtain the following 
recursion relation for the coefficients q 

q 2i -l I 2 
Ci = ^TT 973T ^ °i-mCm, % > 2, Ci = 7 , (17) 
^ L m=l 
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which at q — 0, 7 = 1 generates Bernoulli numbers I?2i, Q = 2 2 *(2 2 * — l)Bzi/(2i)\. One may 
say that (pT7|) defines g-analogs of the Bernoulli numbers [-Bj] ? . Equation (|T7|) works well for 
all values of q. At g < 1 it describes g-deformation of the hyperbolic tangent, since at q = 
one has 

W' + W 2 = <y 2 } W(x) = 7 tanh 7 x, (18) 

which is one-level (soliton) superpotential associated to the Rosen-Morse problem. At q > 1 
one has g-deformation of the trigonometric tangent which is recovered in the limit q — > 00, 

W' - W 2 = 7 2 , W(x) = 7 tan 7 x. (19) 

This superpotential creates an infinite- level Poschl- Teller potential U\{x) with the restricted 
region of coordinate definition: —71 < 2 / yx < tt. On this finite cut Uq(x) = presents an 
infinitely deep potential well. If one sets 7 = simultaneously with q or g _1 then conformal 
superpotentials, W(x) = ±l/x, are emerging. Finally, at q = 1 one gets a standard harmonic 
oscillator problem. 

If q 7^ 0,l,oo, there is no analytical expression for W(x) but some general properties 
of this function may be found along the analysis of Ref. [||], where it was proven that for 
q < 1 superpotential is positive at x = +00. In this case required normalizability condition 
is fulfilled and relation ([TBI) with upper signs really corresponds to physical spectrum. 

At q > 1 the radius of convergence of the series defining W(x) is finite, r c < 00. From 
inequalities 

7 2 _ g 2 — 1 g 21 — 1 

we have < c- < q < c\ , where ci^ are defined by the rule ( |P7|) when g-factor on the 
r.h.s. is replaced by 7 2 /w 2 and 1 respectively (c[ 1,2 ' ) = Ci). As a result, 1 < 2 r yr c /ix < 00/j, 
which means that W(x) is smooth only on some cut at the ends of which it has singularities. 
From the basic relation (|T4]) it follows that there is an infinite number of simple "primary" 
and "secondary" poles. The former ones have residues equal to —1 and their location points 
x m tend to ix{m + l/2)/j, m G Z, at q — » 00. "Secondary" poles are situated at x — 
q n x m , n G Z + , and corresponding residues are defined by some algebraic equations. We are 
thus forced to consider Shrodinger operators ([3]) on a cut [— Xi,Xi] and impose boundary 
conditions ip^(±Xi) = although the potential U (x) is finite at x — ±07. The structure 
of W(x) leads to iJ)q(±xi) = 0, i.e. Vo°' > i s true ground state of H . Note, however, that 
the spectrum E n for such type of problems can not grow faster than n 2 at n — > 00 in 
apparent contradiction with (|16|). This discrepancy is resolved by observation that already 
W\(x) = qW(qx) has singularities on the interval [— xi,x\] so that only H and Hi are 
isospectral in the chain (||). Hence, the positive signs case of ( |16D does not correspond to 
real physical spectrum of the model. 

The number of deformations of a given function is not limited. The crucial property 
preserved by the presented above g-curling is the property of exact solvability of "unde- 
formed" Rosen- Morse, harmonic oscillator, and Poschl- Teller potentials. It is well known 
that potentials at infinitely small and exact zero values of a parameter may obey completely 
different spectra. In our case, deformation with q < 1 converts one-level problem (|I8"D with 
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Eq = — 7 2 /2 into the infinite- level one with exponentially small energy eigenvalues fll6f) . 
Whether one gets exactly solvable potential at q > 1 is an open question but this is quite 
plausible because at q = oo a problem with known spectrum E n = 7 2 (n + l) 2 /2 arises. 

In standard dynamical symmetry approach Hamiltonian of a system is supposed to be 
proportional either to Casimir operator or to polynomial combination of the generators of 



some Lie algebra R IT]. As a result, energy eigenvalues are determined by rational functions 
of quantum numbers. This means that one does not go out of the universal enveloping alge- 
bra. g-Deformation of the universal algebra works with functions (exponentials) of generators 
and, as it was announced, accounts for the presented exponential spectra. 

Indeed, substituting superpotentials (|8J) into relation (|GJ) one finds raising and lowering 
operators 

V& oc H = \{A+A~ - ^4 7 2 ), 

1 — q z 

A + = q 1/2 (p + iW(x))T q , A' = q-^T-^p - iW(x)), (20) 

For real q and 7 the operators A^ are hermitian conjugates of each other. Eq. ([TJ]) insures 
the following g-commutation relations 

A-A + -q 2 A + A~ = 7 2 (l + g 2 ), H A ± = q ±2 A ± H . (21) 

Introduction of formal number operator 



\nH /E 
In q 2 

completes the definition of g-deformed harmonic oscillator algebra in the particular form 



iV = ^f^, N^=n^\ {N,A±}=±A± (22) 
In q z 



The quantum conformal algebra su q (l, 1) is realized as follows [12] , 



[K , K*\ = ±K±, [K + , K-] = - q - q , (23) 

qz _ g z 

i.e. it is a dynamical symmetry algebra of the model. Generators K are parity invariant 
and therefore even and odd wave functions belong to different irreducible representations of 
su q (l, 1). 

In order to generalize basic equation ( |14D we introduce an additional parameter s into 
the superpotential, W = W(x, s), and assume that T q in is a generalized shift operator 

T g W(x,s) =W(qx + a,s + l), (24) 

where q and a are parameters of affine transformation. Although A + is not hermitian 
conjugate of A~ any more, we force them to obey g-oscillator type algebra 

A~A + - q 2 A + A- = C(s), A ± C(s) = C{s ± 1)^, 
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where C is some function of s. Resulting equation for the superpotential 



W'(x, s - 1) + qW'(qx + a, s) + W 2 (x, s - 1) - q 2 W 2 {qx + a,s) = C(s), (25) 

may be called the generalized shape-invariance condition (cf. 0). 
We define a Hamiltonian H as follows 

H = ±A+A~ + F{s), q 2 F{s)-F{s-l) = lC{s), (26) 

where finite-difference equation for the function F(s) is found from the braiding relations 
HA ± = q ±2 A ± H. Now it is easy to generalize formula ([H|). Suppose that a wave function 
ipo, A~ipo = 0, is normalizable. Then a tower of higher states ip n oc (A + ) n ipo gives energy 
spectrum 

n 

E n = F(s) + \ J2 q 2(t ~ l) C(s + i) = q 2n F(s + n), (27) 
i=i 

which can be found by purely algebraic means. If for some n = N normalizability of 
ip n is broken then H has only N discrete levels. In the above presentation we chose the 
simplest form of s-parameter transformation under the action of T g -operator. One can easily 
generalize formula (P7| ) for arbitrary change of variable s in (^), s — > f{s). 

To conclude, in this paper we have described an exactly solvable quantum mechanical 
problem where quantum algebra su q (l, 1) acts on the discrete set of energy eigenstates scaling 
their eigenvalues by the constant factor. In the original version of differential geometric 
applications of quantum Lie algebras an underlying space was taken to be non- commutative 
("quantum plane") and deformation parameter q was measuring deviations from normal 
analysis (see, e.g., Ref.||13||). Here we have commutative space and standard one-dimensional 
quantum mechanics but the potential is very peculiar. It represents g-deformation of exactly 
solvable potentials so that the spectrum remains to be known but it acquires essentially 
functional character. 

It is interesting to know the most general exactly solvable g-deformed potential. One 
of the approaches to this problem consists in repetition of the trick described in Ref . [|Tj| . 
Namely, one can take as particle's wave function a g-hypergeometric function multiplied 
by some weight factor. This would correspond to the transformation of g-hypergeometric 
equation to the form of standard Schrodinger equation for some potential. Another path to 
g-deformation of known models is given by the eq. (|25|) which may have solutions generalizing 
those found by the old factorization technique at q — 1, a — 0. 

Two final remarks are in order. First, affine transformations appearing in fl25|) may be 
used for the definition of g-deformed supersymmetric quantum mechanics JIB] . Second, at 
complex values of g one has meaningful dynamical systems which are exactly solvable when 



q is a root of unity fl6 



The author is indebted to A.Shabat for acquainting with his paper prior to publication 
and for relevant remarks. This work was supported by the NSERC of Canada. 
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